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I Sl'MMAKV  OK  TUKBULI'NT-KI.OW  MODKLS  BASED  ON  THE  REVIEW  BY 
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INTRODl’CTlON 


Sovore  aorodynam  i i-  and  propiilsive  i lUi'rai'l  inns  occur  in  the  fliKi't  of 
VTOl,  aircraft  during  takoolf,  landing,  and  tltc  Lransiti(tn  to  cruise.  In  the 
hover  mode,  unpred  i c I ah  1 e forces  can  he  produced  on  the  airframe  through  the 
interaction  of  the  lift  jets  and  the  gri>und,  a so-called  ground  effect 
problem.  Tlu'Se  forces,  which  occur  as  the  result  of  the  flow  induced  bv  the 
entrainment  of  the  propulsive  lift  jets  and  by  the  de.flected  lift  jets  them- 
selves, c.in  be  either  positive  or  negative.  In  the  former  case  additional 
lilt  is  provided,  but  in  the  latter  case,  a nominally  designed  propulsion 
svsti'm  may  not  have  sufficient  thrust  for  an  adequately  controlled  takeoff. 

The  task  of  determining,  either  experimentally  or  theoretically,  the 
aerodvnamic  and  propulsive  interference  effects  of  the  lift  jet,  the  air- 
fr.ime,  and  the  ground  is  not  a simple  one.  The  numerous  variables  which 
affect  the  flowfield  make  it  economically  difficult  to  study  an  aircraft  con- 
figuration through  wind  tunnel  tests  alone.  On  the  other  hand,  the  flow  is 
so  complex  that  a riuorous  analytical  treatment  of  realistic  cases  is  not 
possible  at  the  present  time. 

fhe  objective  of  the  VTOL  aircraft  flowfield  investigation  in  progress 
at  the  McDonnell  Douglas  Research  Laboratories  (MDRL)  is  to  predict  and 
compute  the  flowfield  assoi-iated  with  a two-dimensional  (planar)  lift  jet  in 
ground  effect.  This  includes  a description  of  the  effect  of  turbulent  jet 
entrainment  on  the  induced  flow  and  on  the  static  pressure  variation  over  the 
u|)per  surface  from  which  the  jet  discharges.  The  two-dimensional,  t ime- 
.iveragid  Navier-Stokes  equations  are  presently  tractable,  and  the  knowledge 
gained  through  their  solution  for  the  planar  lift  jet  in  ground  effect  may 
l.iter  be  applied  to  the  more  complex  three-dimensional  flowfield  problem. 

loward  achieving  this  objective,  MDRL  has  undert.iken  a combined  theoret- 
ical and  exper  iment.a  1 research  effort.  The  initial  analytical  phase,  which 
was  sponsored  by  the  Oft  ice  of  Naval  Research  (ONR)  under  the  subject 
contract,  is  a solution  of  the  two-dimensional,  time-averaged  conservation 
equations,  in  lonjunction  with  a one-e(]uat ion  turbulence  model,  to  describe 
the  turbulent  flowfield  associated  with  a single,  planar,  incompressible 
lift  jet  emanating  from  a static  upi)er  surface  and  impinging  on  the  ground. 

A [>l.'inar  lift  jet  permits  compul.at  ion  of  the  vectored  configuration  since 


1 


Llu'  flowliolcl  romains  Lwo-d  i iirmis  i ona  1 , whereas  llie  ax  i s ymmet.  r i e jet  l)ecomes 
t lirei‘-(.i  imens  i ana  I when  vectorial. 

Tlu’  initial  phase  of  the  MDKh  experimental  effort,  which  was  accom- 
plished under  the  Mcnomu’ll  Donplas  Independent  Research  and  Development 
tlRAD)  rrop,ram,  consists  of  the  accpiisition  of  detailed  llowfield  measure- 
ments lor  the  planar  |i‘l  to  establish  t lu'  validity  c)f  the  computed  flow 

variables.  Altliough  dat.i  have  been  reported  for  a single  axisymmetric  jet 

1-4 

impinging  on  a ground  plane  , an  extensive  set  of  measurements  for  the 
jilanar  configuration  has  not  been  reported. 

This  report  describes  the  analytical  ti’chnique  that  is  used  to  predict 
the  turbulent  flc'wfields  associated  with  the  ])lanar  1 i f t- j e t /a  t r f rame/ground 
i nt I’r.u  t i on . Computed  flow  properties  are  presented  as  functions  of  upper- 
surface  geometry,  Reynolds  number  (based  on  jet  exit  properties),  and  jet 
height  above  ground.  Conclusions  are  drawn  with  regard  to  the  accuracy  of  the 
flowfield  model  and  the  comput a t i ona 1 technit|ue,  and  suggestions  are  made  for 
t ut  lire  St  udv. 
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J.l  Till'  I'li'W  Tniil  i|^uiMtion 

I'lu'  pl.in.ir  impin>;iMK  jft  I lnw  ol  inli  ii'sL  in  tlin  prt'SiTit  .inalyKis  is 
shown  siliom.it  i I'.i  I I y in  Ti^;uio  1 lor  iho  unvooioroil  ooni  i p,ur,ii  i on  . T'ho  joL 
oxits  from  ,i  slot  ol  width  I)  in  i oontourod  nppor  surl.ioo  .i  dist.anoo  II  .ibovc 
tho  p.ronnd  pi. mo.  U'hon  tho  jot  inniiiiKos  norm.il  to  tho  ^;roiind,  tho  roKion  of 
oom|nit  .1 1 i on  I'xtonds  .i  dist.inoo  W on  ono  sido  ol  tho  I i lU'  ol 

svmmotrv.  Whon  tho  jot  is  vootorod,  tho  ri'gion  of  intorost  oxtonds  .a  dis- 
t.inoi'  W on  both  sidos  ol  tho  lino  ol  i;t'omotrir  svmmotrv. 

fho  jot  llowlii'ld  oan  ho  dividod  into  throo  roKiont'.:  .1  froo-ji-L  ro>>ion 
in  whioh  t lu'  I 1 ow  is  ossonti.illy  Lho  s.inio  .as  Lh.aL  of  Jot  issninn  info  an 
nnhoundod  rollon  of  tho  s.imo  fluid;  tho  impinj;omont  ro^ior.,  in  whioh  Lho  flow 
oh.mKos  diroot  ion  with  .1  l.irHO  prossiiro  gr.idiont;  and  tho  w.ill-jot  roRfon  in 
whioh  tho  Mow  tr.ivorsos  tho  surf.ioo  with  zoro  prossuro  Rradiont.  Tho  fluid 
surroundiiiR  tho  jot  is  ontr.iiiu'd  .at  tho  houndarios  of  tho  jot  in  all  throo 
roRions,  o.ausinR  othorwiso  st.it  io  fluiti  to  ho  sot  into  motion.  On  Lho  up\ior 
honndinR  surl.ioo,  this  motion  rosults  in  loo.illy  ri'duood  statio  (irossuros  .itid 
.lorodvnamio  loads. 


d.2  fho  f'.ovorninj^  K([u.it  ions 

Tor  tho  oon  I i Rii  ra  t i on  of  FiRuro  I,  .1  riRorous  fluid  dyn.iniio  .inalysis 
o.innot  ho  m.ido  hv  usiiiR  puroly  i nv  i so  i d- 1 I ow  oaloul.ition  L oohn  i qiios  , ovon 
with  tho  .iddition  of  onipirio.il  or  hoiind.i ry- 1 .ayo r oorrootions.  In  f.iot,  .1 
turhiili'iit  hounda  r V- 1 avi' r proooduro  .ilonp,  tho  Rround  surlaoo  o.innot  simply  ho 
[i.itohod  into  ,111  invisoid  sohonio  sinoo  thoro  aro  prossuro  Rr.idionts  in  hoth 
ooordinato  diroi  t ions  with  no  siiiRlo  dominant  diroot  ion  ol  I low.  Tropor 
tro.itmoiit  ol  till.'  prohli'in  rotpiiros  ;nt  oiliptio  solution  ol  tho  oomploto  I low- 
tiold  with  viscous  ol loots. 

In  till'  prosont  .ippro.ioh,  tho  t imi'-avor.iRod  oontinuitv  .ind  moniontum 
(.N.ivior-.Si  okos)  o<|u.itions  for  stoadv,  t wo-d  i mons  ion.i  I , inoomprossihlo  flow 
aro  usod  to  dosoriho  tlio  moan  motion  ol  tho  fluid.  As  a rosulL  t>f  Lho 
avoraRiiiR  proooduro,  unknown  turhuli'iit  stri'ss  Lorms  ariso  in  thi'  t imo- 
.ivor.iRod  momontum  i'(|u.i  t i ons . To  so  1 vo  lor  tho  turhulont  shear  stress,  a 
t iirhii  I out -k  i not  i o-onorRy  o(|u.ii  ion  is  usi'd  in  oomhination  with  a oonstitutivo 
ofpial  ion  that  rolati-s  tho  sipiaro  root  of  tho  turhulont  kinotio  onorRV  to  tho 
t urhu  I out  V isi’os  i t y . 
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Ill'll  Si- I'v.i  t ion  ol  monii-ntiim: 
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riu-  St  ross  ti-nsor  lor  i noonip  ross  i h U-  flow  is  givon  by: 
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Kqii.it  ions  (2-1)  through  (2-3)  ;irt>  written  in  terms  of  the  primitive 
v.iri.ihli'S  of  veloritv  eomponi-nts  n.,  statie  , ressure  p,  mass  liensitv  p, 
and  moleenlar  viscosity  jj , using  a space-fixed  reference  tl’roiigh  whicli  the 
fluid  flows.  The  overbars  in  the  equations  denote  tli.it  the  qu.intities  are 
d imens i ona 1 . 

To  write  the  t i me-.ive  raged  forms  of  liquations  (2-1)  and  (2-2)  to 
describe  a turbulent  flow,  the  usu.il  Reynolds  decomposition  is  used  i.i 
which  the  i ns  t .in  t .ineous  flow  variables  are  expressed  as  tlu-  sum  of  a n. -an 
component  and  .1  fluctn.it  ing  component. 

"i  ■ <“i>  * "i 

l>  = <(  ^ 

. = <^T. .)  + ■; . . (2-6) 

1 1 ^ 1 K 11 

The  symbol  ^ ^ denotes  a time-averaged  component,  and  the  symbol  ' denotes 
a I 1 uc  t ii.i  t i ng  com[)onent  . 

With  substitution  of  Kcpi.it  ions  (2-4)  through  (2-6)  into  Kquat  ions 
(2-1)  .and  (2-2),  the  following  time-averaged  conservation  ecpnitions  ;ire 
obt  a i ned  : 


(^iiiisi'i  v.i  t i I'll  ot  ni.iss  : 


•<=■> 


Consrrv.K  ion  of  motiu’nt  iim: 


Tlio  forin  of  Lho  l im(J-.ivor;igt.‘d  equation  tor  Llie  eonsirval  ion  ot  mass, 
liquation  (2-7),  is  tlie  same  as  that  of  tlie  instantaneous  equation.  file  form 
ot  the  t i mt'-a  vl*  rajtecl  equation  for  the  i-onservat  ion  ot  momtuitum,  I'.(]uat  ion 
(2-H),  ditti'rs  t rom  the  i nstantaneous  form  tlirouKli  ttie  aiipearanee  of  tlie 
Reynolds  st  ri’ss  tensor,  - p^ujuj^.  .\s  a result  ot  the  latter,  the  system 
ot  f.i]uat  i^'iis  (2-7)  and  (2— h)  is  not  cliisetl,  Theri  arc*  one  sc'alar  c*quation, 
K()nation  (2-7),  and  one  vectc.r  ecjuation,  K(|iiat  ion  (2-H),  hut  there  are  the 
tollowinu  unknowns:  one  scalar, one  vc’Ctor,  1 one*  tensor, 

2.2.2  fu  rhu  1 en  t -K  i neJ_i_cy-Kne  rgy  Kipi.i  t icni 

There  are  five  basic  t urlni  1 eiice-mode  1 i ng  techniques  which  can  be  used 
to  close  the  svstc*ni  of  liquations  (2-7)  and  (2-8).  fhese  are  described  in 
dc-tail  bv  Revnolds  in  a currc*nt  review  article  and  are  summarixed  in 


file  most  lamiliar  rc*presentat  ton  of  the  turbulence  is  the  so-called 
ze ro-f(|ua t ion  mode*!  in  which  the  Reynolds  .strc'.ss  ten.sor  is  refated  to  ,a 
turbulent  (eddy)  viscosity,  which  in  turn  is  related  directly  to  the  time- 
averaged  Velocity  field,  generally  using  a mixing  length.  I'his  approach  has 
had  considerable  success  in  the  calculation  of  boundary  layers,  as  c'videnced 
bv  the  results  of  Cebec i and  Smith  . 

However,  the  xero-ecpiat  ion  model  de  t e r i o ra  t i*s  in  accur.icv  lor  flows 
with  s i g.n  i I i c.int  turning,  high  entering  turbulence  levels,  stagn.it  in  point 
rc*g  ions,  .'ind  bi.und.iry  layer  sepa  r.it  ion  . In  the  problem  under  cons  i de  r.i  t ion  in 


Il)i'  work,  .ill  ol  t lu'so  I low  plu’iiomona  must  bi'  t.ikcii  into  .'iccnunL. 

For  this  ro.ison,  .1  moro  ronipleto  ono-i'(]u.it  ion  modul  o)  t u rbu  1 ciu-t‘  has  boon 
■ iclojitod.  In  subsiMjui'nt  work,  a t wo-taju.i  t i on  model  may  bi'  .adopt  ed  if  the 
sini'lo  i'<|ii.it  ion  proves  to  be  uns.i  t i s I ae  t o rv  . It  should  be  noted,  t IioukIi  , 
that  Reynolds  makes  the  t o I 1 ow  i np,  statement  in  bis  ri'vii'W  article  of 
Relereiice  "However,  it  mav  lu'  that  one  c.an  do  bi'tter  with  Ibi.s  sort  of 

one-e(|u.it  ion  moda'l  in  most  Mows  of  i n t (.■  nss  t , lor  it  m.av  be  easier  to 
specitv  the  1 eiipt  b-sc.i  I e distribution  than  to  comiuite  it  with  a partial 
difliTenti.il  equation.  This  would  be  (la  r t i cii  I a r I y true  if  the  length  scale 
really  should  be  governed  by  the  global  features  of  the  flow  through  .an 
i n t I’g  r.a  1 -d  i t t eren  t i a 1 eipiation.  Hence,  further  stiuly  of  one-equation  models 
is  encour.iged . " 

TABLE  1 SUMMARY  OF  TURBULENT  FLOW  MODELS  BASED  ON  THE  REVIEW  BY  REYNOLDS^ 


k - 


•P.ifiMl  pfTti.jl  r«)u.itiort 

I 

i 


oprr  0)43  2 
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Model  type 

Essential  feature 

Advantaye 

Disadvantage 

1 Zi'to  cquiition 

Uses  only  .1  pde ' for 
the  mean  velocity 
held  and  no  /xle's 
lot  tuibolence 
quantities 

Relatively  simple  and 
therefoie  extiemely 
useful  in  enqineei inr) 
analysis 

Fails  to  handle  some  im[)oitant 
effects,  such  as  strong  surface 
curvature  and  fieestieam  turhu 
fence  and  is  inaccurate  near 
seiiaration  and  in  lioundary 
layers  with  strong  accelerations 

II  One  eun.ition 

In  iKldil'On  !o  a [kJo 
for  the  rnodn  vflocity 
field,  Lisfjs  an  addi 
tional  pde  u?lfM»‘d  to 
tfi«?  tufl)uU*nc<; 
velocity  scale 

Desci  ibes  the  dynam 
ICS  of  the  turbulence 
kinetic  eneiqy  and 
theieby  is  more  accu 
late  than  model 
tyiie  1 

Ret|uiies  knowledge  of  the 
turbulence  length  scale  van 
ation  throughout  the  field  and 
letiuiies  more  com[iuting  time 
tlian  does  model  type  I 

III  T WO 

In  addition  to  a pde 
for  the  meanveiocity 
field,  uses  a pde 
related  to  the  tuif)u 
fence  velocity  scale 
and  a pde  related  to 
the  tur’o'ilence 
lerujth  scale 

Eliminates  the  need 
for  siiecifying  the 
tui bulence  length 
scale  as  a function  of 
IKisition  throughout 
the  flow 

Fails  to  [iroduce  sufficient 
anisotropy  in  the  Reynolds 
stresses  for  general  shear  flows 
and  leriuiies  more  computing 
time  than  do  model  ty(x."S  I 
and  II 

IV  Slf».*ss  fM^utition 

In  addition  to  a (J<le 
lot  the  mean  velocity 
field,  uses  tide's  fot 
all  comiionents  of  the 
tuibolent  stress 
tensor 

Uses  exai  1 equations 
for  the  turbulent 
stress  tensor  which 
are  derived  fiom  the 
Navier  Stokes 
e()uations 

Presents  difficulties  in  satisfying 
all  the  constraiius  which  the 
terms  that  are  used  to  compute 
the  turfiulent  stress  tensor 
must  satisfy 

V LdKjt*  »*(l(ly 
simuldtiofis 

Uses  comimtation  of 
the  tfiree  dimensional 
time  di'iienrlent 
laiqe  eddy  structure 
and  a low  level  model 
for  the  small  scale 
tut  bulence 

E liminates  the 
number  of  rmiiiir  real 
constants  which  are 
le(|uired 

Involves  an  uncertainty  in  the 
effect  of  the  small  scale 
tuifnilence  modeling  on  the 
large  eddy  structuie  and 
retiuirt’S  relatively  large  com 
lulling  limes  for  ()eneral 
shrrai  Hows 

1 ho  ossoiu'i'  ol  t ho  ono-i-qu.it  ion  modi' I is  Lo  dorivo  I rom  tho  insL.'mt.'i- 

noons  'nomoiUnm  oonso  r v.i  t ion  o<|u;Uions  a partial  d i I 1 1- rin  t i a 1 oqnation 

di'Si-r  ih  in>;  tho  Lurhnlont  kinotio  onorov  ol  tho  llow,  k = ( I / 2 ) u ! a ' > . Tho 

^ I j/ 

lattor,  thi'ii,  is  n-latod  through  a oonst  i tut  i vo  oquation  to  tho  tiirhiilont 
visoositv  and,  in  turn,  to  tlu‘  Kovnolds  stri'ss. 

In  tho  pri'Sont  work,  t hi'  t n rbu  I ont -k  i ni' t i c-ono  r>;v  ocpiation  lOrmnlatod 
by  Woltshtoin^  is  nsod.  To  arrivo  at  this  oqnation,  it  is  first  nocoss.ary 
to  dorivo  an  I'qnat  ion  for  tho  kinotio  I'lior^ty  of  tho  moan  motion  bv  adding 
till'  soal.ir  produot  ol  ^n.^  and  tho  t i mo-avo  ragod  momontnm  oqnation  for  ^ 

to  till'  soaiar  prodnot  of  < n.  > and  tho  t i mo-a vi‘ ragod  momontnm  equation  for 

<T>- 


. /r, 


<T> 


- 


9<r>> 


<-‘i> 


"(^''ik)  - K~‘i"'k>) 


, I .,  )>  - p/uin’  )) 


(2-9) 


Whin  Kq'.iat  ion  (2-9)  is  oontraotod  (i  = j),  tho  result  is  tho  k i no  t i o-ener  gy 
oipiation  of  tho  moan  motion: 

_ 9 d/2)<(ri  )><(n  )><(ri  )> 

„ ^ i = - < n,>  - 


ik>  - 


(2-10) 

)<'n.'> 

\ \i/  /-  \ \i/ 


i 


I 


■ 


To  oht.iin  I lu'  t‘(|u.U  ion  lor  the  kiiu'tio  oiutkv  oI  tlu-  turbulent  lluctua- 
tions,  the  senl.ir  prodmt  ot  u.  .ind  t lu'  momentum  ecjiKition  for  u.  | Kquat  ion 
(d-J)l  is  atliied  ti'  till’  seal. If  produet  of  u.  .ind  the  moment  urn  equation  for  u. 
[Kc|uation  (d-2)].  Kqu.it  ions  (2-  t)  through  (2-b)  .are  introduced  into  the 
resulting  etjuation,  whiih  is  t ime-a  ver.iged , combined  with  Kqu.ation  (2-9), 
and  ci'ntracted  ti'  give  t lu“  mean-kinetic-energy  etpiat  ion  of  the  turbulent 
fluctuations  (the  so-called  t urbu 1 ent-k i ne t i c-energy  equation). 


(2-11) 


Uv  rewriting  various  derivatives  and  hv  introducing  the  conservation 
of  mass  equ.ition  and  the  defining  equation  for  K<|uat  ion  (2-11)  can  he 

shown  to  be  equivalent  to  the  form  of  the  t urbu 1 ent-k i ne t ic-ene rgy  equation 
usi'd  by  Wolfshtein^: 


The  terms  in  Kqu.ation  (2-12)  have  the  following  physical  meanings: 


Dx 


Convection  of  turbulent  liinetic 
energy  by  the  me.an  motion 


9 


I’rodiict  ion  o|  t iirbii  1 fn  L kincti( 
energy 


l urhulfiU  din  us  ion  oi 
LurhiilunL  kinetic  enerj’v 

V'iscoiis  dilfusion  oi  turbu- 
lent kinetic  energy 


Viscous  dissipation 


2.2.  ) Wo  I t slit  (^i  n ' s One-Kfj^uat  ion  Turbulence  Model 

Id  ordi'r  to  establish  a relation  between  the  various  turbulence  proper- 

t ii's  and  the  turbulent  kinetic  energy,  the  turbulence  mocieling  procedure  of 

Wolfshteii/  will  be  followtul.  It  is  assumeil  that  a scalar  turbulent  vis- 

1 os  i t V u , c.in  !u'  tie!  ined  bv 
t urb 


- K''!";) 


‘‘t  urb 


wlu‘ri‘  the  turbulent  viscosity  is  computed  I roni 


i|  , 

I urb 


l/2_ 

c tk  V . 
h M 


(2-li) 


(2-1 '4) 


In  Ktiuat  ion  (2-\^*)  c is  a <-onstant  (hdermined  I rom  i-xper  iinent  , c 
J p 11 

and  9 is  the  length  scali*  lor  viscosity. 

U 


0.22. 


10 


Ill  .uKlition,  it  is  .issumi'd  tli.it  the  tiirhiili-ni  diltusioii  ot  tiirbiilenL 
kinetie  t'lu  ri’.y  e.m  he  represi-nt  ed  liv 


y",",  > 


9x.  \ k,iurh  dx 

) ' ) 


Will'll  i'.qiiat  ion  (2-lS)  is  .idded  to  tlie  vi.scoiis  ilil  fusion  of  turbulent  kinetic 
- 2-  - 2 

energy  ti(()  k)/(Dx.‘'),  tlie  result  is 


_3  / ^^turb  iik_  \ “ 5^k 


3x . \ (7,  , 3x . 

j \ k,turb  j 


3x . 9x . 


k , t iirb 


By  del  inin>>  an  effective  transfer  coefficient  for  the  turbulent-kinetic- 
energy  diffusion. 


_ ]j_  ^ turb 

k.eff  ” a,  C 

k k,turb 


Kquat ion  (2-Ib)  can  be  rewritten  as 


-turb  I ok 
)i  + I 

'^k,  turb  / 9x  . 


k , e f f 


where  O,  is  unity  and  0,  , is  a constant  determined  from  data,  0,  ^ , 

k k.turb  k,turb 


Final  Iv,  it  is  assunu'd  that  the  visca)us  dissipation  can  be  represented 


by 


i 


I 


is  a I'oiisLant  evaluated  from  data,  = 0.42,  and  is  the  length  scale 
for  d iss  ipat  ii'n. 

When  K(|uat  ions  (2-13),  (2-18),  and  (2-19)  are  substituted  into  liquation 
(2-12),  the  result  is 


’‘turb 


9x . 
J 


<=j> 

9x . 

,1 


+ 


J) 

9x. 

J 


V 


I) 


(2-20) 


K<|uat  ion  (2-20),  the  so-<’alled  Wolfshtein  one-equation  turbulence  model, 
is  written  in  terms  of  tlu-  mean  velocilv  field,  turbulent  kinetic  energy, 
and  expi- r imenta  I 1 y determined  constants  and  length  si’ales.  The  latter 
will  be  specified  when  particular  geometries  are  ctuisidered  in  the  analysis. 

2.2.4  Non-Oimensionjil  Forny_(2_f  the  Time-Averaged  ContinuUy,  Momentumj_and 

lurbiilenl -Kinet  ic-i;neri;y  ITpiat  ii'ns 

I'.qu.'itions  (2-7)  atid  (2-H)  are  written  in  x,  y coordinates  with  the 
associ.ited  u,  v vi- 1 oc  i t V i'omi)onen  t s . I'he  lontinuity  ecpi.it  i('n  li.is  been 
used  in  writing  the  momentum  equations  in  the  non-conservation  form  shown 
be  low. 


i 


C('iiKiTv;it  ii'n  of  mass: 


- ()_  (2-21 

9x  9y 

CiMiscrvat.  ion  of  monuaUiim  I L Lmo- a vo ra>;o(J  Navier-:'tuk.es  or  Reynol  ds  o(]iiaL  i ons)  : 
x-momont urn  oquaLion: 


j ay  9v  9x  3y 


y-raomoiUiim  oqiiat  ion: 


9y  9y 


A 


Di't  iniii>;  I'ciiiat  ions  lor  tlio  turbuloiit  and  ollo.  livo  viscosities: 


>f  r 


|i  + u 


t nrli 


(2-24) 


nrh 


11  ;i 


(2-2'i) 


Wo  1 I'slu  e i n ' s one-equation  turbulence  model,  Kquat'on  (^-20',  takes  the 
tol  lowing;  form  in  the  x-v  coordinate  svslem; 


With  t lie  specification  of  the  length  scali'S,  the  syste-m  ol  equations  is 
closed,  subject  to  t lu‘  rec|uired  boundary  conditions. 

Kquations  (2-21)  through  (2-2b)  are  normalized  by  introdiuing  the 
following  dimensionless  variables; 


(2-27) 
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(2-28) 


V 


o 


■f  f 


'efj 

f)V  D 
o 


(2-29) 


u 


t iirb 


turb 


pV  I) 
o 


(2-30) 


Tbi?  dimtTisionli'.ss  t'fju.U  ions  are  ^iven  below. 
Conservation  of  mass: 
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(2-31) 


y = 


(2-32) 


i. . 

(l/2)pV^ 


(2-33) 


k = - - 

- 2 
V 

o 


(2-34) 


^ = 0 . 
t)x  3y 


(2-35) 


x-:nomi-nt  um: 


u 


V 


HI 

iV 


' ‘P  + 

2 )x  'efl 


'Hi  , 'Hi 

ju  JJii'H  III rl) 

Jv  ''iy  Hx  ()y 


V-  momi'iU  um : 


♦V 

’V 

U * + V 

♦ V 

'■i’  4-  u 

'V  ft  f 


411 

(ju  Lj.i^r^l') 

(jy  dx 


+ 


1 , , 

IV  1 urb 
• x 4X 


•Hi  , 
-)v tiirb 

Hy  TV 


WolfslUuiu  oiu'-uquat  ion  turbulence  model: 


d k d k 

+ V „■  = , . , 

• lx  IV  turb 


+ 


-IX 


,ef  r 


Ik 

IX 


) 


+ 


>v 


k.eff 


'V 


) 


(2-36) 


(2- 37) 


3/2 
' 1) 

(2-38) 
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V'lU’  ft' 


U 


cif 


I 

( urn 


f 2-  VI) 


k.ut  t 


U,  i<v 

k 


^“,urb 
k , t uri) 


(2-40) 


I V I) 

Ku  = . (2-4  J) 


2.2.')  Voj^t_i  I'J.  t_N’  / St  ruam-Fuiu.'  l iun^  Koriri_  o I lu*  T i mi-'-Avfr.'igotl  (k)iu  i mi_i  tv  , 

Momunt  um,  aml_  rurj)u  1 lmiLjK  iju' t u;j-Knc K(j_uatJons 
Kriu.nt  ions  (2-)'))  llirougli  (2-27)  aro  solvod  in  terms  of  t i mo-ave  r.iged 
vortii'itv  and  stream  (unction,  which  are  defined  (ly  (die  t'ol  Inwini;  vquutions: 


dx  iiy 


(2-42) 


I'i 

ay 


Ll 


(2-43) 


Ox 


V 


(2-44) 


Tile  relations  lietween  tin'  dimensionless  ami  dimension.il  variables  (i)  and  IJ) 
•are  given  by 


(V  /O) 
o 


I 1 


(2-4S) 


V 1> 

n 


Thf  vi'rtii'iLv  iiMiisporl  c(|u.uion  is  licr  i vod  l>v  d i t t ori'nt  i;iL  iii^  liquation 
(2-\h)  with  rcspi'ct  to  v , di  I t oront  iat  in^  li(|uat  ion  (2-17)  w i t li  ruspoct  to  x, 
and  subtraoting  tiu'  tornior  I roir  tlu-  latter. 

Vortieitv  transport  efpiation: 


( 1 + Re*  ) 


Re 


O 


+ Ki'*  ) 


9 


+ Re 

av“ 


•)y 


1 

9 

•) 

')  ■} 

■} 

2 

2 

9 

7 

. 1 

o. 

r- 

Ul 

(2-47) 
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X f)  y 

9x3y 

^ ';-2 

9 + 9 

3x“  dy 

-> 

;)y " 

> 

■)x‘’ 

0 

:jx"' 

7 

^y^  - 

Ihe  tol  lowing  substitution  has  bi-en  niatle  in  tlu'  |>revious  equation  lor 
s imp  1 i e i t y ol  not  <i  t i on  : 


►- 


*\urb 


(2-48) 


Conservat  ion  ot  mass  is  ensuri'd  through  ilid  inition  of  the  stream 
function,  which  is  di'termined  froiii  a Poisson  i-quat  ion  as  obtained  by 
combining  l.quat  ions  (2-42)  through  (2-44). 

Poisson  I'qiiation  lor  stream  1 unction: 


iT.'k  = 

2 

by" 


(2-49) 


I 8 

t 


I 


P!quation  (J-iH)  is  ,ilsi>  wrilti'ii  in  ti  rms  nt  ..  ,,  ,in<l  • . 


I 


i 


? 


f 


I 

I 


Wo  1 1 sht  i' i n t iitiuil  onoo-moiir  1 miu.it  inn  : 


(2-50) 

Tlio  boundary  conditions  on  (j,  and  k wliicli  arc  applii-d  to  Equations 
(2-47),  (2-49),  and  (2-50)  arc  spocitiod  wlu-n  s|H‘<  itic  >;e('ine trios  are  con- 
si  do  rod . 

Onco  ttio  vorticitv,  stroain  fum-tion,  and  turbulont  viscosity  variations 
aro  computed,  tl\o  velocity  component  distributions  and  the  static  pressure 
li  i s t r i but  i ons  must  be  ca  1 1.  u i ati'd  . Tbe  velocity  components  are  readily  ob- 
taiiu'd  t'rom  tlu’  detinition  of  the  stri'am  lunctit'ii.  Equations  (2-4'j)  and 
(2-44).  Ihe  pressure  I ield  is  computed  throuKh  a solution  of  the  Poisson 
equation  for  i>ri-ssiire.  This  e(|iiation  is  di'rived  by  differentiating  Equation 
(2-!b)  with  respect  ti>  .k,  tl  i f ferent  iat  ing  Etjuat  ion  (2-17)  witlt  respect  to  y, 
and  adding  the  two  ecjuations. 

Poisson  I'quation  for  static  [)ressuri': 


(2-51) 
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K(|u.it  ion  is  siilvnl  lollowiin',  I hr  soliilion  ol  K(|  ii.i  L i ons  (2--'<7),  (2-V) ) , 

.iiui  (2-')0);  tli.U  is,  l.lu'  ri^;ht"li'"''l“-‘^  ion  is  ;i  known  cjuniUity. 

riu'  voiiij' 1 ox  i 1 V nssiu  i .1 1 oil  willi  tlio  solution  ol  tho  I’ois.son  o(|u.ilion  lor  stalir 
prossuro  is  t li.U  llio  ma^nitudo  of  t ln'  prossnr<'  is  not  known  al  tlio  boundaries. 
In  till'  prt  si'iU  anaivsis,  the  prossuro  1 'ovi' 1 is  sot  at  tlii'  luid  |)oint  ot  -a 
lioundarv,  and  tlu’  iinniHuit  urn  o(|u.it  ion  doscriliinn  tlio  pn-ssuro  Kr.idient  along 
till'  boundarv,  I'itlu'r  Ikiuation  (2-22)  or  F.quation  (2-2'J),  is  integrated  to 
obtain  tile  tioundarv  stat  ii’  pressure  vat'iation.  I'be  oonstraints  wliioh  are 
imposed  oTi  tlU'  remaining  liound.ir  i I'S  aiu'  t lu'  normal  |>ressure  gradients  given 
bv  till'  two  raomi’iUum  equ.it  ions. 


i. 


I hi:  m'mkkicai,  sdi.iition  schkmk 


Hk’  I'Dupl.il,  clliplir,  I'.i  r t i .1 1 -il  i r t'cTcMi  t i .1 1 t i ons  that  de.scribf 

thi'  ll<'V  .iro  non ! i lUM  r , and  a nnmo  r i la  I scIutho  is  lu'Cfssarv  to  olitain  <i 
solution.  To  si' I vt*  (he  tl<n»’field  efjn.itions  for  a ji‘t  discdiarging  from 
.1  onrved  uppi’r  surface,  a eonlormal  ma|)pini’,  techniciue  is  used  to  map 
till-  irreKular  phvsical  plane  into  a rer  tangu  I ar  eomput  a t iona  I pl.ine.  Tlie 
>;ovi’rn  in,n  eipiaticms  are  then  rewritten  in  terms  of  tlie  computat  i<jna  1 plane 
roordinati'S  and  si' 1 ved  usinR  f i n i t e-d  i f f erence  procedures. 

In  this  section,  tlie  conformal  mapi'ing  and  finite-difference  techniques 
are  described. 

1.1  The  Conformal  Mapping 

lo  solve  the  governing  equations  lor  a I low  witli  a contoured  upper 
boundary,  whicb  simulates  the  lower  surface  of  a fuselage,  an  inv<  rse 
conformal  mappinv,  procedure  is  introduced  in  the  analysis.  Details  ot  the 
mapping  proceduri'  and  its  application  to  the  flowfield  equations  ;ire  given. 

1.1.1  helper  ipt  ion  oj’  the  Mapping  Procedure 

flu-  mapi'ing  technicpie,  which  was  originally  devised  at  MDRI.  by  Holfman, 
is  best  I'xplaiiU'd  with  reference  to  Figure  2.  Initial  Iv,  a f in  i t e-d  i f t e rence 
comput.it  ion.i  I plane  with  coordinates  (',0)  specitied.  The  distance 
between  nodes  in  the  ■■  direction  is  a and  in  the  q direction,  b,  where  a and 
b .ire  not  neiessa  r i 1 y e(|u.il.  Stretching  functions  are  then  introduced  in 
I'.ich  coordinate  diriction. 

.i=f,(-)  (3-1) 

‘ = !,,(•)  . (3-2) 


With  thi'se  rel.it  ions,  .i  mapping  plane  (lJ,A)  is  determined  which  permits  liner 
resolution  of  the  flowfield  in  regions  where  t lu>  gradients  ol  the  computed 
vari.ibles  are  severe.  Fin.illv,  a conformal  m.apping  given  by 
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Figure  2 Inverse  conformal  mapping  with  stretching 


riu'  bomul.irv  i-oml  i t ions  imposi'd  on  llriii.i  t i mis  (5-4)  ;imi  (5-'))  to  I low  I rom 
(ihvs  u'.i  1 const  r.i  i nt  s wiuti  t lu’v  .n  i‘  known  ,it  t ho  boundaries  and  t nim  integra- 
tion ot  till'  I .nub  V- 1\  i i-mann  rol.it  ions. 


d _ cb'  -jy 

d..  '■  d>  iVi 


( 3-b) 


d',  bx  _ _ ^ 
d‘  o d .i  d-', 


(3-7) 


wlu’ii  tlio  X .and  t lie  y bonnd.irv  distributions  are  not  known. 

riio  unif|uo  toaturo  of  tbo  prosont  mapping  srliomo  is  tiiat  instead  of 
spooifving  llio  i-oonl  in.at  i>  li  ist  r ibut  ions  in  the  pliysical  plane  and  accepting 
wh.ati'ver  computalion.il  pl.ine  results,  ;i  numerically  convenient  computational 
pl.iiu'  is  specilied,  .ind  the  corresponding  computational  mesh  within  the 
specilied  physical  bound.iries  is  computed.  for  this  reason,  the  mapping 
scheme  is  .in  inverse  procedure. 

Kqu.it  ions  (3-4)  .ind  (3-3)  are  discretized  in  unexpanded  form  using  the 
ci-nt  r.i  1 -d  i 1 I ereiicc  , finite-difference  .approximations  of  the  derivatives  given 
in  Appendix  A .and  the  point-ot -the-comp.ass  grid  notation  defined  in  figure  3. 
llii-  following  lorms  ot  K.ipl.u-e's  etpi.ition  result: 
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whert-  the  coot  t i c i elll  s .ire  lefilU'd  by 
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Finite  difference  grid  parameters 

b Moimali^ed  grid  dimensions  in 

f and  r)  directions  respectively 
a I)  Ratio  of  grid  dimensions 
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Figure  3 Five  point  finite-difference  stencil 


Kill'll  Immul.irv  ilislrihut  ions  ol  x .mil  v nro  nlu.iim-il  I rom  I ho  (!.m  hy- 
Kii'inann  ro  1 .1 1 i ons  , I'.tiu.il  ions  ( i-(i)  ,iiul/or  { )- / ) nro  i n t oj',r.  1 1 oil  tor  x ainl/or  v 
using  a soiotui-ordo  r-aoo  n ra  t o (|uailralnro  solionio.  Ki|u.it  ion  ( !-8)  is  so  I vi-d 
1 irst  for  t ho  v field  using,  point  relaxation,  and  I lion  h(|nation  (i-4)  is 
solvod  lor  till'  X tiold  using  point  ri' luxation. 

In  ordi'r  to  iK'torniino  t ho  aoonr.iov  of  the  nuiiiorioal  iiia])ping  proiodnro, 
a ti'st  oaso  without  ooorainate  strotohing  was  tried  lor  whioh  the  relation 
botwi'on  the  mapp  i ng- p 1 aiu'  and  phys  i oa  I -p  I ane  ooordinatos  is  known  oxaotly. 

Iho  prohloin  lormulation  and  the  funotional  ri'lnlions  In'twoon  (x,v)  and 
are  givi-n  in  Kignro  -i(a). 

In  the  nunii'rioal  solution  ol  l.aplaoo's  equation  lor  y,  y is  s])ooiliod 
analvtioallv  along  the  upper  and  right  houndaries.  fhe  rem.iining  boundary 
distributions  of  v and  all  the  boundarv  d i s t r i bn t i ons  of  x are  determined 
from  pbvsieal  eonst  mints  or  the  tiaur  bv- K i emann  relations.  The  latter  are 
used  instt'ad  of  the  anaivtie  boundarv  values  in  order  to  di'termine  the 
.u-euraev  of  the  rpiadratnre  sebeme  in  eomputing  the  boundary  eoiulitions. 

figure  4(b)  shows  the  eomiiarison  between  tin.'  exact  and  numerical  mapping 
for  the  test  case  with  ;■  = ' / i . With  a 41  x 41  f i n i t e-d  i I t ereiice  grid,  the 

analvtical  and  numerical  ma|>pings  agree  at  worst  within  10 


I 

i 


Figure  4 Conformal  mapping  test  case 
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3.1.2  App  1 i c.U  ion  of  t ho_  Mpppi  n^J’roceduro  i o i he  J’pvt'rn  inj^^  J;^uat  i on.s 

Ihe  t.  ime-.iviT.i^i'd  oonservaL  ii'n  e<|viat  ions  ami  I urb\i  1 em  e-niode  1 equ.ation 
derived  in  .SeiLion  2 which  deserihe  a two-d  imens  i I'lia  I i nccimpress  i h 1 e flow  can 
he  written  in  tlie  Poisson-  or  t ranspor  t-eciiiat  ion  forms  given  below. 

Poissc'n  ecination: 


•)  T 

n + ill:  = 

0 ^ 9 

;3y‘' 


(3-15) 


Fi'r  the  Poisson  equation  for  stream  function,  liquation  (2-49), 


= ilJ  (3-16) 


J,,  = - w , (3-17) 

and  for  the  Poisson  eqnat  icni  for  static  [iressure.  Equation  (2-51), 


P = P 


(3-18) 
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(3-19) 


Transport  eqnat ion; 
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(3-20) 
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Ki'r  thi-  tr.inspi'rt  oquation  lor  vortioity,  Equation  (2-47), 


= 10 


T j = 1 + Kc‘0 


')  ■)  0 9 

Sy  Sx”"  ')x  9y^ 


and  lor  t ho  t ii  rhii  1 onco-modo  1 oquation. 
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tainat  ion  (2-’0), 
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Applving  t hf  fontormal  m.ippin^;  proiH'dure  with  coord  i n.i  te  stroLching 

dose  r i lu-d  in  Section  d.I.l,  derivatives  ol  a general  variable  '}j  in  the 

phvsical  plane  (x,v)  can  he  expressisl  in  terms  of  di>riva  fives  in  t hi’  corapn- 

tational  pl.uie  (F,.r\)  with  the  formulas  for  f D'{))  / (3x)  , (3‘''r)  / (3x“)  , (d'Ji)/(3y), 

2 2’ 

0"'|')/(')y  ),  and  ( 3“  f ) / ( 3x3y ) given  in  Appendix  K.  Using  these  representations 
of  the  derivatives  in  conjunction  with  the  (lane  hv-R  iemann  relations,  Kr|uations 
(3-1'j)  and  ( i-20)  assume  the  following  torms  in  the  comput a t i ona  1 plane; 
I’oisson  equation: 


In  these  equations,  1}  is  the  mapping  modulus  and  is  delined  by 


,i : 


ii^y 


n-ri) 


Q 


Squat  ion  (j-!3)  is  not  i-onvenicnl  for  tiio  mimfrical  calculation  of 

2 

Q~  because  (x,y)  are  solved  in  terms  of  (^,,ri).  Ifliat  is  needed  is  Q in 
terms  of  the  derivatives  |(9x)/(;iC)<  (0x)/(9q)|  or  ii'v)/(a^J,  (9y)/(9’i)l 

which  can  be  directly  written  in  finite-difference  form.  i'lirough  a[)pli- 
cation  of  the  chain  rule  and  the  Cauchy-R i emann  relations,  it  can  be  shown 
that 


.r  = ; — r-~2  • 

/ii  . hii  '*211 

ydii  9r  j + \^dA  h^| 


Similarly,  the  mapping  derivatives  (D;i)/(9x)  and  (9A;/(3x)  appearing 
in  Ixiuat  ion  ( 3-  32)  can  be  expressed  in  the  forms 


9 a ^ 2 d r 

9x  ' du 


(2-35) 


and 


9 A 2 dq  i)x 
9x  ” dA  9p 


( 3-  36) 


2 " 2 2 

As  shown  in  Appendix  15,  the  di'r  i vat  i ves  (9  0 ) /(9x“)  , f i"f)/(9y  ),  and 

(9^1 ) / (9x9  y)  which  appe.ir  in  tlie  right-h.iinl  sides  of  Kquat  ions  (3-Jl)  and 

2 ,2  .2  .2 

(3-32)  reqviire  the  calculation  of  (9  ;i)/(9x‘')  and  (9  A)/(9x  ).  The  latter 
ari'  computed  with  the  following  formul.as  using  numerical  differentiation: 
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( 5- 58) 


Kin.illv,  I lu‘  far  l I's  i an  vnlncilv  cnmpoiU'nt  s arc  cominiLea  from  Llio 
Tf 1 at  ions 
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5.2  'Hu'  [•' i n i t f-l)  i I I ortMui’  liThnirpu^ 

Ttu'  i’oisson  ('(Illations  for  stream  tiinclion  ,uul  static  pressure  e.in  be 
solved  mime  r i ca  I 1 V witliont  aitlieiiltv  using  tlie  conventional  eentral- 
aitterence  algoriiln.  This  is  not  the  case  for  the  vorticity  transpon. 
eipiation  and  the  Woltshlein  1 iirhii  I ence-inode  1 transport  e(|uation.  The 
coetticieiits  ol  the  1 irst-order  derivatives  in  t hesi  equations  contain  the 
Kevnolds  number  as  a mu  1 t i p I i ca t i ve  factor,  and,  as  a result,  with  the 
standard  cent  fa  1 -d  i M ereiice  algorithm,  the  discretixed  system  of  ecpiations 
is  diagonally  dominant  tor  onlv  a I imited  rangt'  in  the  magnitudes  of  the 
(•Of  t f i c i en  t s . Diagonal  dominance  is  necessary  to  obtain  (unvergence  in 
the  iterative  solution  of  the  discretized  system  of  ecpiations.  Dt'tails  of 
the  t i n i t e-d  i t I e rence  analysis  and  tin'  iterative  i rocedure  for  solviirg  the 
disiretized  ecpiations  are  given. 

5.2.  1 Disc  ret  izat  ion  ot  the  (loverninj^  Kcpiat  ions 

The  Poisson  ecpiation,  licpiation  ( 5-51),  is  discretized  using  the  standard 
cent  ra  1 -d  i f t erence  .ip[)  rox  i m.i  t i ons  given  in  Appendix  A and  the  point-of-the 
compass  notation  detined  in  l■'igure  5. 
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I'll  i ssmi  o(|u.i  t ion  : 
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whoro  till'  I’oo  f f i i' i on  t.  s aro  ovaluatod  I rom  I ho  lol  lowing  rolaLions: 
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Prior  to  I ho  d i so ro t i /a t i on  ot  t ho  transport  ocpiat  ions  lor  vortioity 
and  stream  tnnotion,  Ihpialion  (i-!2)  is  rowrillon  in  the  torm 
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( !-49) 
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In  t lio  ptvsi'iit  work,  Kf)iiation  ( i-4H)  is  writton  in  t' i n i l o-d  i t‘ I oronco 
form  iisinK  tlu’  oii-ci'tit  r.i  I -d  i f foronoc  (A(d))  nigorithm  dovolopi'd  bv 

Hoffman  at  MDKI.  . I'hf  cssimui-  t>f  this  miHliod  can  ho  illiistratod  by  consider- 
ing till'  derivative  of  (i^,‘)/(3c)  of  Ikpiation  (1-48).  Using  the  fivo-pitinl 
diftoronco  stinoil  shown  in  I'ignro  ) and  point-of-tho  ct'iiipass  notation,  this 
derivative  can  he  ovaliiati'd  at  point  I’  using  the  following  truncated  Taylor- 
serii's  represiuit  at  ion  and  standard  ciuit  ra  1 -d  i t f <■  riuico  .ippri'x  i ma  t i on  to  the 
f i rst  dor  ivat i ve : 
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In  the  AUl)  scheme. 

t he  der ivat  i ve  ( d 

is  r. 

terms  of  lower-order  derivatives  hv  d i I 1 1' riit  t i a t i ng 
re.spei  t to  f).  The  derivative  ( ''*')/(  ’ i)  in  l.f|u.i'  ion 
an  analogous  f.ashion  with  the  AUl)  algorithm. 


( i-4  5) 


o i nod  ind  is  prossi'd  i n 

K(|ua  t i on  ( i-48  ) with 
(’.-48)  is  ropii'seutod  in 
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With  this  fitiiti'-ai  M‘c-rc'iu-f  niftlKul,  Kqu;ition  (J-48)  has  the  cl  i sc  re  1 1 zfd 
tc>rm  yi'on  below, 
iranspoi  t ecpiat  ion: 


(i-54) 


ll 

+ 9A 

T 6 ^ 

\ ’ 

p he 

P IT 

p :)t| 

P P dr/r)r, 

( ^ 

:i.s 

+ 

4 A 

V 9' 

P 

P hri 

p p 3f,,-)n 

(3-61  ) 


2 


C3-62) 


Boiiml.'irv  (.'oru]  i t.  i ons  ari_‘  i lU  rocliicod  into  K<|ii;i  1.  i ons  (3-41)  and  (3-54)  in 
till'  I'ol  lowing  minnor;  WIkmi  a Dirichlft  oondiLion  is  sper  i f iod , Llu>  known 
v.ilui'S  oi  ; art'  .'mposod  on  iho  <1  i sc  ro  I i zod  ofjuat  ions  ovaliiated  at  the  grid 
lino  adjaoont  to  t ho  boundary  of  intorost.  When  a Noumann  condition  is  im- 
posod,  till'  normal  dorivativo  is  d i sc  rot  i zod , and  t ho  values  of  ^ at  the  grid 
lino  adjacent  to  the  boundary  outside  the  solution  domain  are  related  to  the 
values  of  ; at  the  grid  lino  adjacent  to  the  boundary  inside  the  solution 
dom.iin.  rho  discretized  o(|uations  evaluated  at  the  boundary  points  are 
included  in  the  system  of  equations  to  he  solved. 

Till'  voloi'itv  components  are  computed  I rom  Hquat  ions  ( 3-39)  and  (i-40) 
using  conventional  cent ra 1 -d i f f e rence  finite-difference  approximations. 


5.2.2  Solution  ' tj^yj)  i sc  re  t i zi'd  Kqua  t i I'lis 

Ihe  1 i n i t e-d  i 1 I erence  eipiat  ions  are  solved  iti'ratively  using  point 
relax.it  ion  with  the  caliulating  si'(|uence  illustrated  in  Figure  5.  It  all 
the  comput  .It  i ons  were  performed  with  .a  .single  computer  code,  the  required 
I'omputer  storage  would  be  excessive.  Consequently,  three  sep, irate  codes  .ire 
useil  in  series. 

With  the  first  computer  code,  the  m.ipping  I'lpiat  ions  .iri'  solved  intera- 
t i ve  1 v for  V and  x using  ;i  ri'l.ixat  ion  f.ictor  of  l.b  until  a loc.il  maximum 
residual  between  successive  iter.it  ions  of  10  is  .achieved  in  the  field.  The 
m.ap[)ing  coordin.ites  .ind  derivative's  ari’  stored  for  use  by  the  second  computer 
I'ode  . 

With  the  second  computer  code,  the  flowfield  equ.it  ions  .ire  solved  itera- 
tively for  the  stri'.'im  function,  turbulent  kimtic  energy,  and  vorticity  using 

respective  rel.axation  f.ictors  of  l.t),  0.8,  and  O.H  until  maximum  local 

-4  -')  -4 

residuals  betwei'ii  successive  iterat  ions  of  lo  , 10  , and  10  , respi-c  t i ve  1 y , 
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.Ilf  .iihii'Vi'cl.  (It  .1  l.iinin.ir  .snlui  ion  is  licsiroii,  ilic  tiirhulent  viscosity 
is  si't  to  zero  evi‘ rvwiieri’  in  t lie  lielil,  .iiul  t lu-  solution  ol  the  tiirbulont- 
k i lU' t i e-i'iie  r^v  isjii.ition  is  hypasseii  in  the  iti’ialion  loop.)  I'tie  mapping; 
coord  i n. It  es  and  di-rivatives  .iiul  the  stre.im  lunetion,  turhuUMit  kinetic 
i‘ner);v.  ,uul  vortieil-.  il  i s t r i but  ions  are  stored  for  use  hy  the  third  computer 
code  . 

With  till'  third  lomputer  codi‘,  the  Poisson  equat  ion  lor  static  pressure 

is  solved  iti'fativelv  usin^  a relaxation  factor  of  l.b  until  the  maximum 

-4 

residual  in  tlu'  field  between  successive  iterations  is  10  . In  the  final 

set  ol  comput  .1 1 ions , the  velocity  com[ionents  are  cominited. 

All  computer  codes  <ire  written  in  I'ORTRAN  IV  for  use  on  the  Control 
Data  Corporation  CYKKR  173  computer  of  the  McDonnell  Douglas  Automation 
Comp.iny,  St.  Louis.  , 


Cdlciil.itK)f<  of  Ww.  ro.ippjfu} 
cooidinati?s  atul  (l«*r ivdtiv«?s 


iHK  coMi’iiTKi)  ri.oWF I i;u)  soi.rrinN's 
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In  t lu‘  prusinit  work,  tlowriold  solutions  h.ivo  boon  carriod  out  for 
thri'o  Ki’>'mi’trio  variations  iif  tlio  planar  impinK>nK  jot  shown  in  Figuro  1. 

I'ho  first  goomotrv  is  a t ro<‘-uppi'r-boundary  oon  f i gura  t ion , so  oallod  booauso 
tho  tippor  bound, irv  of  tlio  solution  domain  is  sovor.il  diamotors  bo‘ 1 ow  the  slot 
from  which  tho  jot  d i sch.argos  . This  coni  i gura  t i on  wms  soloctod  since  it  h.is 
boon  .in.ilvzod  by  Wolfshtoin^  and  offers  .a  mo, ins  of  comp,iring  the  MDRL 
solutions.  Tlu'  second  geometry  is  a par.i  1 I o 1 -p  1 at  e cont  i gur.it  ion  in  which 
the  jot  disch.irgos  from  ,i  slot  in  tho  upper  surlace  ,ind  impinges  perpendicu- 
l.ir  to  t lu‘  lower  surf.ico.  This  con  f i gur.i  t ion  was  included  in  the  studv  since 
it  otters  ,i  mi'.ins  of  evaluating  tho  tlowtioid  model  ,ind  solution  scheme' 
without  till-  ,idded  numi'ric.il  com|)lexity  of  conform.il  m.ipping.  i'he  third 
geometry  is  the  con  f i gur,i  t ion  of  prim.iry  interest  in  which  the  jet  issues 
from  .1  I'urvi'd  plate  th.it  simul.iti’s  .i  tuselage  undersur  f.ice . Both  .in  un- 
viHtored  and  .i  viutored  jet  are  considered,  .and  in  the  formi-r  c.ise  solutions 
.ire  presented  for  two  j et -he igh t -.above-ground  siuicings  .ami  for  two  Revuolds 
numbers  with  oiu-  of  the  11  v.ilues. 

. 1 file  Free-L'pper-Boundary  (leomet  ry 

flu‘  f ree-uj>per-bound.ary  geomi’try  is  i 1 1 us  t r.a  t I'd  in  Figure  b.  A tree, 
turbuU'Ut,  symmetric  plane  jet  issue.s  vertically  downw.ard  onto  tb.e  ground 
[il.ine.  Since  the  problem  is  geomi't  rica  1 ly  symmi’tric  with  respect  to  the 
centi'rline  of  t hi'  jet,  only  h.a  I I of  the  flow  dom.ain  nei'd  be  solved.  The 
vorticity  .ind  stre.am  fumtion  .are  antisymmetric  .about  the  centerline,  .and 
the  turbulent  kinetic  energy  is  symmi'tric.  Ihus,  the  solution  dom.ain  extends 
.1  dist.imi'  11  .above  the  ground  pl.ane  .and  .i  dist.ance  W to  the  right  ol  the  line 
ot  geometric  symmetry.  .No  conlorm.il  m.ipping  is  retpiired  since  the  boundaries 
are  rect.angul.ir,  and  coordinate  stretching  is  not  used.  Thesi'  s imp  1 f icat  ions 
result  in  x = f,  and  y = '!• 

I he  imposed  bound. iry  conditions  are  shown  in  Figure  b.  On  the  upper 
boundary  ol  the  solution  dom.ain  within  the  ji't,  the  imposed  velocity  profile 
determines  the  stre.am  fuiict  ion  .and  vorticity  distributions,  and  the  turbulent- 
k i ne t i I -ene rgy  profile  is  specified.  Outside  the  ji't,  along  the  remainder  ot 
the  upper  boundarv,  I lu'  flow  is  taken  to  lu'  irrol.it  ion.al  .ind  normal  to  thi' 
free  surface  with  no  turbulent  fluctuations.  fhe  boundary  conditions  on 
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Figure  6 Free-upper-boundary  geometry 
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I lu‘  i«'t  i-.-n  1 1 r I i 111'  rt’sii  1 I Irom  svmnu'irv;  on  I lie  lower  ;;iirlac-o  tlicy  result 
1 rom  I lie  no-slii>.  imperme.ih  1 e wall  eonslraiul.  'I'lie  f I owr  at  t lie  riKlil 

houmlarv  is  t akiui  to  exit  (i.arallel  to  the  wall  w 1 I li  no  ;;  racl  i lui  t s of  turbulent 
kiiu'tie  energv  in  tiu'  mainstream  ilireel  ion.  Tlie  dimensionless  lenr.th  scale 
variations  used  are  those  |)ro|>osed  by  Wo  I I sli  t e 1 n ^ . 


i - exp  ^-Re.\^^k'^‘-'v  j 


(4-1  ; 


- ex])  ^ -KeA^^k  ' j 


(4-2) 


whiTe  A and  A,,  are  empirical  Iv  determiiK'd  constants  with  values  of  O.Olb 
U 1) 

and  0.21,  fi'spec  t i VI' 1 V . I'ollowlng  Wollshtein,  it’  tlu'  magnitudes  of  tlie 
length  si-ales  excei'd  0.1,  f and  art.'  set  to  0.1. 

The  I 1 ow  varialiles  were  computed  with  II  = I,  W = I,  and  a Reynolds 
numlier  of  1)  POf)  based  on  jet  properties  at  the  up|H‘r  boundary  of  t he 
solution  ilonia  i n (denoted  "i"  slati“  in  Figure  b)  , which  corresponds  to  a 
Kevnolds  numlier  of  II  000  based  on  properties  at  tlie  noxxle  exit  plane 
(denoted  "n”  state  in  i-'igure  b).  Computations  were  perlormed  using  a 
■'ll  X ■'*  I unilorm  I 1 n 1 te-d  i 1 1 eri-nce  grid. 

Contour  plots  of  tlie  normalized  primary  flow  variables  are  given  in 
Fig.iiri'  7.  rile  vorticitv  il  I s t r i bu  t I on  shows  the  convection  ot  ui  to  tile  riglit 
and  the  development  ol  tlie  lower  wall  houmlarv  laver,  and  tlie  stream 
limction  plot  illustrates  t tie  I'li  t ra  i nmi'Ut  of  fluid  into  the  I ree  jet.  Tlie 
distribution  of  t ii  rbu  I on  t - k i no  t ic-energy,  which  iias  a peak  value  ol  0.04 
imposi'd  at  till'  u|ipi'r  lioundarv  within  the  ji‘t,  sliows  t lu‘  convection  of  k 
toward  tile  rig, lit  bound.iry  .ind  a [ilateau  of  maximum  k in  the  vicinity  of  tlie 
stagnation  point,  indiialing  a rate  ol  liigb  tiirlnilence  generation  in  tliis 


■'lO 


reg,  I on  . 


Tlu’  norma  1 i zt“ii  primitivo  I I ow  v.ir  i ah  I os  art*  shown  in  Figure-  H.  Tht* 
contour  plot  ot'  the  x component  ot  velocity  clearly  shows  the  development  of 
the  wall  houiularv  layer  aiul  the  entrainment  of  tluid  into  the  jet,  and  the 
plot  ot  the  V component  ot  velocity  depicts  the  decay  in  v as  the  stagnation 
point  is  .iiipro.iched . The  variation  ol  static  pressiiri*  shows  a small  normal 
pressure  gradient  throu>;h  the  wall  houndary  layet  and  a general Iv  tavorahle 
pressvire  gradient  along  tht*  wal  1 as  tht*  flow  acceleratt*s  tt^  tht*  right  . 

The  1 ree-upper-bountiarv  geomt*try  defined  in  Kigure  h has  been  an.ilyzed 
bv  Kolfshtein^  for  the  cast*  of  laminar  flow  with  tht*  results  [)resenti*d  in 
Reference  9.  A comparison  of  the  MDRl.  laminar-flow  soluti('ns  with  those  of 
Kolfshtein  for  11  = W = 1 and  Reynolds  numln*rs  of  1 and  1000  is  given  in 
figure  9.  Wol f shtc in ' s t urhu 1 en t- f 1 ow  solutions  are  not  for  a solution 
domain  with  11  = W = 1,  so  the  corresponding  comparison  for  the  cast*  t>f 
turbulent  flow  is  not  applicable. 
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I’l.  p.i  r.i  1 1 1 ■ 1 -p  1 . 1 1 1'  p.ccnu  t rv  i - shown  in  lip.nri-  lit.  A i n rhn  I oni  , 
svniTulrii'  pl.u\>.-  ji't  d i ‘".Un  rpcs  t rom  i slot  in  t ho  njUKT  ll.it  pl.it  i-  .iinl 
impinpis  iii'fm.il  to  t ho  Kfound  jil.ino  .i  dist.inoo  i holow.  A-.  in  t ho  Iroo- 
iippo  i-hotind.i  r I'.isi  , t ho  t 1 ow  is  so  1 vod  on  I v lor  t ho  h.i  I ! [>l  ino  (oxlondiiiK 
.1  disi.inoi  l\  to  I ho  rivtht  ol  tin.'  ji't  ooitt  or  1 i lU’ ) to  i.iko  .idv.in  t .im-  of 
■•'.ooiTio  I r i o s vinmo  t r . '■  ■ nin,  no  oonlortn.il  m,ip|)inn  is  noodod  sinco  i ho  honnd.i- 

rios  .iro  roo  t .inp.u  1 nr , and  ooordin.ito  strotohiiiK  is  not  nsod;  thus,  x = 

•md  V = • . 

rho  hound. irv  oonditions  lot  this  m.'omotrv  ,iro  illiisir.itod  in  l•'iKuro 
lO(.i).  Iho  v.'iluos  ot  v'or  t i o i t , stro.im  lunotion,  .ind  turhuloiit  kiiu’tio 
i'iior>;v  aro  spooit'iod  within  t ho  onti-rinn  ji't,  iiul  alon^  tin-  iippor  w.ill  t hoy 
rosult  from  the  no-slip,  i mpo  rmisih  1 o w.ill  const  r.i  i nt  . I'ho  oonditions 
imposoil  I'U  tin.'  romaininK  throo  ho'  nd.irios  ,iri'  idoni  io.il  to  thosi'  iisisi  for 
t roo-upiio r-hound.i r V ^ooniot  rv  . 

For  t ho  p.iral  lol-|)|alo  oonl  iyural  ion,  tho  .loouraov  ol  t ho  .issumpt  ion 
of  no  j;r.idionts  in  tho  1 1 ow  proportios  .it  tho  riyht  hound.irv  do  t o r i o r.i  t I's 
as  tho  r.it  io  H/W  >f  tho  solution  domain  is  iiu  fisisod.  .\  ntoro  >;otioral 
r i yht -bounda rv  oondition  is  oiu'  in  which  a norma  1 -prossuro-grad i oiu  protilo 
and  .1  t tirhu  I on  t -k  i lu-t  i c-oiio  r;.’,v  protilo  .iro  iiii])ost-d.  iho  oorrosponditiK 
vortiiitv  and  st  roam  lunotion  d i s t r i hut  i ons  must  t Inui  hi'  oottipiilod  i t o r.i  t i vo  I 
I rom  tho  I’oisson  orpiat  ion  for  stro.im  funot  ion  and  tho  v-momontum  oqiiation, 
wluro  (ap)/(iy)  is  known.  Ilowovor,  prosontlv  tho  roipi  i ri'd  oxpi' r i monLa  1 
prossu  ro-it  rad  i on  t and  t urhii  I out -k  i no  t io-oiurnv  d i s t r i Ini  t i ons  aro  not  known 
lor  till'  oonl  igiir.it  ion  of  intorost,  .and  tho  .assumption  ot  p.iral  lol  flow  is 
i nvokod . 

For  tho  pa  r.i  I 1 1' 1 -p  I .It  o gi'oiiu'trv,  aotu.il  I luigl  h-so.a  I o distributions 

• iro  simil.irlv  unknown.  Consoipiont  I v , tho  to]  lowing  lanist  r.a  i nt  s hasod  on 

simplo  ordor-of-m.agn  i t udo  .argumonts  h.ivo  boon  usod  to  suggissl  ]ilivsioallv 

ro.i  I i s t i o d i s t r i hii t i ons  o I V,  = .at  . i w.i  1 I , ■ 0 ; at  v = 11/ d , 

U 1) 

I = 11/2  and  (dU)/(dy)  = 0;  and  within  tho  rogion  ol  tho  ontoring  jot  givon 
by  (0  1 X 0.2'j,  v = H),  • = 0.5.  ihoso  oonditions  aro  satisi  iod  hv  tho 

fill  lowing  [Ml  1 ynom  ia  I s : 


I = 9 

I 


(II  F I ) 


(1  + ''ill) 


(.'i-n 


4 ) 


0 _ X L 0.25 


I 


fn  till'  ri’^iiin  ().2'>  ■ x O.'i,  t lu'  follnwiiiK  romh  i n.it  ion  ol  ,in<l  f ^ 
is  iisi'il  to  |■^■pr^■Sl■nl  > : 


= 4 


( 0 . S 


X)  ‘I,  + (X  - O.;”'!)  9. 


0.25  • X 0.5. 

(4-5- 


I'qii.itions  (A-i)  throiiKli  (A-5)  .iro  ontirolv  K‘-‘'’nu't  r i r tuni’tions  .-incl , unlike 
liqu.it  ii'iis  (4-1)  ami  (4-2),  do  not  lont.iin  exponential  terms  that  depend  on 
the  tiirbiili'iu  kini'tie  ener^v,  Revnolds  niitnber,  and  empirical  constants. 

I'iitiire  10(h)  illustrates  the  1 enqt  h-sca  1 e distributions  for  a par.i  1 1 e 1 -p  lat  e 
qeometrv  with  II  = 2. 

(ia  1 cu  1 at  i ons  were  made  with  11  = 2,  W = 4 . bH , atid  a Reynolds  number  of 
10  000  based  on  propertit's  .it  the  jet  exit  slot. 

Ki^ure  II  shows  the  computed  iirimarv  flow  variables.  Vorticity  is 
convected  to  the  ri^ht  boundary  of  the  solution  domain  with  development  of 
,111  .ittaclii-d  bounil.'iry  layer  on  tlu*  lower  surf.ice  and  a separ.ited  boundary 
l.ivi-r  on  the  upper  surface.  (Sep.ir.it  ion  oicurs  on  a solid  boundary  where 
. = 0,  .md  in  the  resent  conf  i>;urat  ion  this  occurs  very  near  the  slot  edge.) 
The  region  ot  rec  i riu  ating  flow  is  cle.arly  illustrated  by  the  stream 
function  plot.  rill'  m.i.ximiim  v.ilue  of  turbulent  kinetic  energy  imposed  on  the 
entering  |et  is  0.04,  ,ind  .i.-..  the  contour  [ilot  of  k reveals,  there  is  a 
✓ dec.iv  I rom  this  value  throughout  tlu'  fiolii.  Over  the  right  half  of  the 

solution  domain,  the  turbulent  kinetic  energy  is  fairly  uniform. 

Ilie  primitivi'  t low  vari.ililes  .are  g,iven  in  Figure  12.  The  contour 
(ilots  ol  the  two  velocitv  components  detine  in  more  det.iil  the  region  of 
1 — rec  i rcii  1 .1 1 i on  between  the  pl.ites.  I he  distribution  of  the  x component  of 

Velocity  illustr.ites  tin’  upper-  and  lower-wall  boundary  layers  and  the 
. eiii  r.i  Inineiit  of  fluid  into  the  jet.  Thi'  sharp  |)ri>ssure  gradient  in  the 

imiiingemeiit  region  is  shown  in  the  plot  of  the  static  pressure,  where  the 
value  of  pressure  at  the  st.ignuion  point  li.is  been  si't  e(|u.il  to  one. 
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Figure  11  Primary  flow  vanaljles  for  the  parallel  plate  geometry 
(H  = 2.  W ^ 4.68,  Re  10  000) 
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Figure  12  Primitive  flow  varidbles  for  the  pHrdllel  plate  geometry 
iH  2,  W ^ 4 (>8,  Re  10  000) 
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Till'  i-ont  orm.i  I m.ipp  i ii)’,  si-lirnu'  di'Sc  r i li«  i|  in  Si'.  I inn  i.l,  willinut  t lie 
nso  oi  s t ret  ill  i np.  Inint  inns,  w.i-.  .ipplicd  in  i Ii  i s p.c’nnntrv,  .iiui  snlntinns 
wnfi'  I'.irrind  nnt  tnr  t Im  I I nw  rinid  i • inns  spnn  i I ind  in  I ivpirn  1 

In  ill  iMSns,  till'  mapping  ntpiatinns  wcrn  snlvcd  Idr  nnlv  t Im  li.i  I I - 
[ilaiii';  in  trnat  tim  nnmplnln  llnwtinld  Inr  I lu'  vimlnrnd  ji'l,  t Im  rnfpiircd 
m.ipping  do  r i va  t i vi‘S  wnrn  rcfiniind  .ilioiil  I Im  I i im  nl  >;cnnm  I r i n svmnintry. 

With  rnitard  tn  t Im  halt  ,1  am,  W dnnnlns  t Im  width  nf  t Im  cnmpu  I a I i nn.i  I 
plane,  .ind  W'  dniintns  t Im  ha  I I -width  nt'  the  upper  hnmid.iry  nf  the  phvsieal 
plane.  The  hall-width  nf  the  mnde I , nnrmalizeii  hy  tim  jet  slnt  width,  is 
•t.hH,  hut  a ennv'erpi'nt  snliit  inn  nl  the  maiipinp,  erpiat  inns  enultl  nnt  he 
aefiii-ved  when  W'  was  extended  tn  this  v.ilim.  The  mappinp,  pmeedure  retpiires 
that  hndv  slnpes  he  finite,  hut  fnr  the  e u r ved-|i  I a t e (’eniiietry,  (dy)/(dxl  is 
extreinelv  large  in  the  vieinitv  nf  x = a . h8 . The  largest  value  nf  W'  tnr 
whieh  a ennvi' rgeii t snlutinn  cnuld  he  nhlained  tieereased  as  II  was  rediieed . 

I'nr  this  reasnn,  the  values-  nl  W'  speeilied  in  Figure  ii  are  nnt  etpiivalent 
tn  the  hall-width  nt  the  eurved-piale  gmimi'lry. 

Tlu'  hnundary  innditinns  in  the  eniiiiiu  t a t i nn.i  I pi. me  impnseil  nn  the 
finwfield  snlutinn  fnr  the  eu  rved-p  1 .1 1 e gentnetrv  are  siinwn  in  Figiin*  14. 

tains ider  It  st  the  hnundarv  enndilinns  Inr  the  unveetnred  jet.  Figure 
14(h).  Onlv  the  right  h.ilt  nf  the  llnwfield  is  snlvral  sinei'  t lu'  prnblein  is 
.ig.'iin  gennel  r ie.i  I 1 v svmmelrie  with  respiu-t  tn  t lu’  jet  centerline.  The 
vnrtieitv,  stream  funetinn,  and  t u rhu 1 en t -k i net i e-ene rgy  prnliles  are 
specified  within  the  entering  jet.  On  the  upiier  and  Inwi'r  surfaces,  the 
hniind.irv  cnnditinns  result  I rnm  the  ennstraint  nf  nn-slin,  imperiiieah  1 1' 
w.i  1 I . On  the  jet  eeiiterline  the  svnimet  rv  ennditinns  .iri'  impnsed,  and  at 
the  right  hnundarv  t In-  I I nw  prnpert  ii'S  iri-  t .iki'ii  tn  he  unifnrni  with  res|ieet 
tn  f.  It  is  recngn  i XI  li  th.it  this  eniulit  inn  heennu'S  i ne  riMs  i ng  1 y in.ieeuratt' 
with  deere.ising  W,  and  in  future  wnrk  exper  inuml  .i  I I v nht.iiiU'd  f 1 nw  (print  i t i I's 
can  he  imposed. 

(iniisider  new  the  tiound.irv  cnnditinns  used  with  the  vectored  jet. 

Figure  1.4(c).  In  this  cnni  igpirat  inn  the  snlutinn  must  he  c.irried  nut  fnr 
the  entire  llnwtii'ld  with  a width  .!W  in  the  cnniput.ilinii.il  plane. 

Within  the  enrering.  jet,  the  vnrtieitv  .ind  stream  liinctinn  v. iri.it  inns 
-ire  cnmpiiti'd  I rnm  the  veincitv  profile  nht. lined  hv  rnt.it  ing  the  prntile 
imjinsed  in  the  unveetnred  case  thrniigh  .111  .ingle  with  respect  tn  t lu' 
line  nf  g,enmetric  symmetry.  file  t iirliu  I ent -k  i net  i c-eiieigv  d i ;.  t r i hu  t i nn  is 
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■ ilso  oht.iiiu'cl  thr!'UKli  .in  .ixis  nU.ilion  nl  l hr  nn  vci- 1 o red  k il  i s t r i hii  I i on  . 
V.ilue.s  of  ,,  Ji,  ,md  k on  I lie  solid  snrl.irt-s  ;ire  ev.iln.iled  1 rom  the  no-s  1 i ji , 
ini|n>rmeahle  wall  eondition.  The  eonsir. lints  ^’,iven  hv  ( i ) / ( if  ) = Qkl/fa") 

= 0 art.'  imjiosed  on  both  I lu'  ri^tlit  and  leli  houmlaries  of  the  solut  ion  domaiti 
however,  the  condition  (;)iJt)/(o-)  = o is  applied  on  1 v at  the  left  houndarv, 
ami  a t^t  distribution  (editained  through  i n t I'p,  r.i  t i on  ot  ,i  speeilied  u profile) 
is  imposed  on  the  rig, hi  hound. irv.  This  is  necessarv  to  I’valuati"  the  lower- 
wall  stream  function  since  ill  can  no  longer  he  i n t i-gr.i  t t'd  along  the  jet 
center]  ini’  from  the  upper  houtidary  to  the  lowi'r  boundary,  as  is  done  i ti 
the  unvectored  case. 

The  hmgth-scale  distributions  used  are  those  given  hv  Kquations  (4-T) 
through  (4-'))  with  11  repl.iced  by  coordinate  of  the  upper  surface. 

Kor  the  unvectored  con f i gu ra t i on , 41  grid  points  were  used  in  both  t he 
and  ’ directions  in  the  f in  i te-d  i 1 fereiici-  solution.  I'or  tin-  vectored 

configuration,  41  grid  points  were  used  in  the  T.  direction  ami  21  in  the  n 

d i rec  t ion . 

4.1.1  IJiTl 

I’rimarv  and-  primitive  flow  variahU's  computed  lor  the  curved-plate 
geometry  with  11  = 4,  W = 4.b8,  and  Re  = 10  000  are  shown  in  Figures  IT  and 
16,  respectively,  for  the  c.ise  of  normal  ji't  i mp  i ngi'men  t . The  flow  profiles 

are  qualitatively  simil.ir  to  those  presentml  lor  the  para  1 1 e 1 -p  1 a t e geometry 

However,  in  t tu‘  c urvml-p  1 a t e con  I i gu  ra  t i on  , with  the  larger  .aspect  ratio  ol 
the  solution  domain  (ll/W),  the  ri'gion  ol  rec  i rcu  I a t i ng  flow  iloi'S  not  lie 
compli-ti'ly  within  the  conlininv,  hound.i  r i I's  . 

Solutions  Wert'  computed  for  the  coni  ii'.iiration  with  Re  = 100  000,  hut 
they  are  indistinguishable  1 rom  those  computed  for  Re  = 10  000.  The  reason 
is  ap|i.irent  1 rom  Kcpiation  (2-14);  for  large  Reynolds  numbers,  the  norm.ilized 
mo  I ecu  I. a r visctisitv  (1/Re)  is  a very  small  contribution  to  the  elfective 
viscosity.  The  latter  is  domin.ated  by  the  turbulent  contribution,  which  is 
controlled  by  the  speeilied  geomi'trii'  length  sc.ales  and  not  by  the  Reynolds 
number  based  on  the  molecular  viscosity. 
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Normalized  x component  of  velocity 


Normalized  y component  of  velocity 


Normalized  static  pressure 


Figure  16  Primitive  flow  variables  for  the  curved  plate  geometry 
(H  4.  W 4 68,  Re  - 10  000) 
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Figure  18  Primitive  flow  variables  for  the  curved  plate  geometry 
(H  2,  W = 3.68,  Re  1000) 
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Fiyure  19  Primary  flow  variables  for  the  curved  plate  geometry 
(H  = 2,  W 3 68,  Re  10  000} 
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Fiqurp  20  Primitive  flow  variables  for  the  curved  plate  geometry 
(H  - 2.  W - 3 68.  Re  10  000) 
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I)  i 1 I I ■ I ■ '111  1'  - hilwi'fii  tin  I'lmpiiliil  Ili'W  |)  ri'|iiT  t i >■  s im  1 lu-  r.i-.c  willi 
Ki-  - llKIO  .mil  Ki-  III  <10(1  ifi  I'l- 1 .1 1 i v>- 1 V sin.i  I I tnr  l lu-  rc.isnn  r i t iil 

pr>  V i>ni>  1 . . lh>  1 ■ v>  1 >'l  l <i  i i>ii  I > ;u  i i ,i  t i i m ni'.ir  I !ic  t .ip.iiM  t i nn  puinl 

i 111  I .'.mis  willi  Kcvh"  I il '■  imiTil'i  t . i-;  i I I us  I r 1 1 I'd  b'.’  l hu  pints  nt  1 u rbu  I nut 

liiiu'lii  niiiri’v.  UnWfvni  , ihiTn  IS  tin  s i )',n  i t ii.mt  •.Mri.itinii  in  l lin  nturain- 

tii-nt  InVi-l  tnr  ,1  i b.mm  in  K<-  i rnrii  Ititii)  tn  1<)  i)i|i). 
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. . . llic  1 1 1 ri'il  II  I 

As  .m  inil  ial  ti'Sl  c .isi*  lur  l lia  cnmpiit  i-r  nulo  ust-d  ti>  coinpii  1 1'  t lie  vec- 
j,i  1 luwt  i i- 1 il , .1  si'lutiiiii  was  rarrii-il  out  lor  an  nn  vor  t oroil  contiH'ira- 
lion  to  vi  rilv  that  in  ilio  limit  iiiK  oasi-  oi  ■■  = f)”  tho  i-orrict  ilow  patti-rn  is 
ot.tainoii.  I'iKiiro  22  prosi-nts  tho  primarv  Ilow  v.iriahlcs  lor  tlio  oomploLi' 

, nrvoil-plato  poomot  rv  with  H = . V - i.hH,  and  Ko  = M)  000.  Siiu'o  tho 

valuo  ot  tin-  1 owor-surt  aoo  siroain  tniution  c;in  in'  sot  to  .'oro  in  this  case, 
t hi  prohlom  was  solvod  subioct  to  N'oumann  honndarv  oonditions  on  all  throo 
prim.irv  flow  variahlos  at  both  tho  loft  and  riKlit  hound. irios.  Th.it  is.  tor 
tlio  nnvo.  torod  jot  tho  rij’ht  bound.irv  oondition  = spooiliod  in 

liipiro  l'.(.  ) w.is  roplaood  hv  ('.l/fh-)  = 0,  .md  tho  low.-r  hound, iry  condition 

- 1 (0)  was  roplacod  l)V  . = 0.  l-ip,uro  22  olo.irly  il  lust  rates  tho 

1') 

.isvminoi  rv  ol  tho  vorticitv  .and  stro.im  function  il  i s t r i hu  t i ons  .ind  the  symmetry 
.i|  the  t urhu 1 out -k i not  ic-onornv  distribution.  The  dittoroncos  hotwoon  tho 
curved-p  l.ii  o tlowiiolds  computed  tor  tlio  lull  pl.ano  llit;uro  22;  .and  tho  h.alt 
pl.inc  (iiciiro  17)  .iro  attributed  to  tho  dittoromo-.  in  tho  i i n i to-d  i I i oronco 
■:/  sh  si.'cs  usi'd.  The  tormor  wore  i alcul.itod  with  .i  co. irsi  r mesh  to  conserve 
mil  hi  no  computing  time. 

A solution  w.is  .ilso  carried  out  for  tin  comploti  curvod-plato  goomotry 
witii  11  - 2,  l>  = = II)  , and  Ko  = ID  dDU.  In  this  vectored  jot  noomotrv, 

■ill-  .'.iluo  ol  tho  stri’.im  function  on  the  Iowa  r surlaco  is  not  knowat  a print  i 
mu  ss  .1  sLiiiplo  horizoiiL.il  moi.iontum  hai.uico  i .ipplii'd  ,issunun>t  invisi  Ld 
impiiu'omont  tlow.  In  lieu  ol  this  .issiimption,  the  houndarv  condition  imposed 
:l  tlio  ri,;ht  boundary  w.is  ov.ilii.iLod  1 rom  the  x-componont -o  f-vo  1 oc  i t y computed 
lor  tile  unvoctori'd  jot  with  tho  s.imo  velocity  prol  ilo.  Iho  ri'sultin^  solu- 
tion w.is  .m  unro.ilistic  flow  p.attorn,  which  indicates  th.it  physio. illy  correct 
■ .'lid  it  ions  must  ho  imposed  .it  the  llowliold  houndar  i o.s . Uork  is  continuinK 
I'll  the  soliii  ion  ol  tliis  I lowl  iold. 


1 
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Sl'MMAKY 


Starting  with  t lu'  i-om|i  i i' I c flliptic  .una.  ■ r v.i  I i . in  iqu.itinns  tor 
i ni'nmpri'ss  i li  1 1> , sti'.ul".  viscnus  tlnw,  an  ana  1 v t i ca  1 tlnwtii-ld  model  lor 
a sinpde  t wo-d  imi'iis  iinia  I (planar)  lilt  jel  in  pronnd  ettCet  has  Inen 
tormiilated  and  solved  mime  r i ea  I 1 v . The  Kovernin.:  erpiations  ire  t hi'  Reynolds 
I'qiiat  ions  aiul  a l u rbii  1 inici'  model  equal  ion  whiidi  are  e.ist  in  terms  o|  two 
transport  ecpiations  (one  for  vorticitv  and  oiu'  1 Or  turbulent  kinetic  eiierqv) 
anil  two  Poisson  etpiat  ions  (one  for  stream  tunct  ion  and  one  tor  static 
pressure).  The  latter  are  solved  usiiiK  I i n i t e-d  i I 1 e rence  procedur.s  in 
ciiniunct  ion  with  a contormal  mappint;  technique  to  treat  solution  domains 
with  niin— I'l'c  t an>tu  1 a r uppia"  houndarii's.  flowt  ields  havi*  been  conijiutA'd  tor 
threi'  geometries  with  various  jet  impingement  con t i gu rat  i ons  and  Kevnolds 
numbers . 

5 . 1 fauic  IjjS  i Otis 

With  regard  to  the  numerical  schemi',  no  significant  problems  havi-  heiui 
encountered  in  solving  the  coupled  system  of  equations  describing  the  time- 
averaged  tiowtield  for  the  unvectored  lift  jet.  For  t hi'  Vi'ctored  I'ontigura- 
tioii,  a better  definition  of  the  r i gh  t -hound.i  r v conditions  is  requ  i ri'd  to 
obtain  a phvsicallv  ri'.ilistic  repri'Siiit  at  ion  of  tlu’  flowtield.  (Convergent 
solut  ions  ot  the  contorm.il  mapiiin'g  equat  ions  without  the  usi'  ol  coordin.ite 
stretching  have  fiei-n  obtained  for  a contoured  u|i|ier  surl.icv  e.xcept  whiui  th.it 
surt.ice  cont.iins  ,i  region  ot  verv  large  sloiie. 

With  regard  to  the  analvtical  model,  the  com|iuli“d  solutions  display  all 
the  expi'cted  cha  rac  t e r i st  i cs  ot  a jet  i mii  i ngi’ment  llowfield;  tin*  lower- 
surtace  boundarv  l.iver  with  the  typical  wall-jet  velocity  protile,  tlu'  upper- 
surtace  bound. irv  layer  with  sep.i  r.i  t i on  , the  considerable  pri'ssure  grailient.s 
and  high  level  ot  turbulent  t 1 uc t uat i ons  in  the  impingement  region,  and  the 
distribution  ot  velocity  along  the  jel  centi'rline.  These  I 1 ow  p.itterns  are 
iiased  on  (lostulated  liound.irv  conditions  at  the  jet  entrance  and  exit  iilanes 
and  on  [los  t u 1 .i  t eil  I engt  h-sca  1 1’  distribiit  ions.  The.se  weri'  deriveil  using 
orde  r-o  f -magn  i t iide  arguments  obt. lined  I rom  sinqile  physical  riMSoning,  since 
the  experimental  dat.i  required  for  a more  precise  detinition  were  not 
.1  va  i 1 at)  I e . 
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llii'  .u  i-ur.u’v  ol  till'  .issiinu’i!  hnuiui.irv  coiiil  i t inns  .md  UsijMli-.si’nli- 
li  i si  r ihiit  inns  will  In-  i-s  t , il>  1 i slu-d  lliroiu;li  ,i  ci'mp.i  r i son  nt  l hi-  t hc-nn-t  i r.i  1 
t li'wl  it-lds  w i t ri  J.it.i  I'lu.iini'd  in  .in  Ml'KI,  IKAl)  u-si  pr«i)’,r.im.  I'lic  l.ittnr 
is  (.ifvnu-d  to  t ho  misisnroraonl  ol  pl.in.ir  )i  I imp  i np,o!iu.Mil  t'lowlioUls  in  iloso 
groumi  tltoit  for  the  I l.it-pl.ito  and  curvod-pl  .it  o goomi- 1 r i o.:  con.s  i do  red 
in  the  present  <-i>ntriit  work. 

'i . ’ Kooomme_nd.i_t  iojis 

Following  a dot. li  led  ooinp.i  r i son  ol  t lu-  oompulod  .and  mo.'isurod  I low 
proin-rt  ii“' . the  hound. irv  .ind  1 ongt  li-so.i  1 o distributions  should  In-  iniprovi-d 
.IS  lu-ooss.i  r-.-  to  provide  .i  mori-  prooisi-  noiiel  ol  llu’  I low.  It  ;i  simple 
par.inietrie  studv  e.innot  In-  u.sed  to  pr>>vide  improved  1 1'liiU  h-sea  1 e v'.iri;it  ion.s, 
.1  t wi'-e<pia  t i on  model  ot  tiirhulenee  is  proposed  to  eompule  nu're  .leeur.ile 
leiigtli  Sc, lies.  The  latter  Would  theii  he  iisi-d  in  the  ont--efpii  t i c'li -mode  1 , 

e 1 in  i na  I i ng  the  .iddition.il  n.iehine  c-c'mpul  .it  ion  lime  recpiired  to  perlorm  the 
e. I leu  1 .It  ions  with  -i  fw.i-ecpi.it  ion  tormul.it  ion  ol  the  tiirhulenee  tield.  I'sim; 
the.'  improved  lengtii  sc.ile-;,  the  tlowtield  solutions  ior  the  i.-ont  i gura  t ion 
ol  primarv  interi'St,  the  eu  rved-p  I .i  t e geonu-t  r-.  , should  he  c-arried  out  tor 
r.inges  in  v, lines  ol  Re  .md  II  in  order  to  est.ihlish  the  ini  I ueiiee  ot  t he 
l.atler  .'II  tin-  tiirhuient  eiit  r.i  i nnient  .ind  the  surl.iee  st.it  ie  iiri-ssure  distri- 
hul  i oils  . 

'rt'ith  till-  completion  ot  this  cll.-rl  lor  the-  c-.ise  ol  i lU'omp  ress  i h 1 e flow, 
it  is  proposi.sl  th.it  the  .in.ilvsis  he-  c'Xtended  to  the  e.ise  ol  e omp  re  s s i h 1 e 
flow  tor  the  s.ime  eon  t i -,',11  r.  1 1 i ons  . Such  .111  extension  is  justitied  since 
in  re.ilitv  the  lluid  in  the  lift  jet  is  at  .1  t empe.  ,it  ure  much  higlier  than 
th.it  ot  the  surrounding  lluid,  .ind  the  imniedi.ile  llowlic-ld  h.is  l.irge  densitv 
r.id  i i-n  t s . For  .1  Vrol.  .liiir.ilt  in  ground  eltc.-ct,  this  e.in  he  signilie.int 
.IS  .1  result  ol  one  Of  more  ot  the  lol  lowing,  me.li.in  i sms : 

(1)  Dells  i I V. iri.it  ions  helwc'cn  tin-  amhieiit  .air  and  the  less  dense  lilt 
lilt  jet  h.ive  .111  ini  luc-nei'  on  the  eii  t r.i  i nment  ot  air  at  the  t rc.'e 
hound.it  ies  o|  the  tree  jet  .ind  the  ground  w.i  I I jets. 

(2)  .Mixing  ot  the  .imhieiii  eiitr. lined  air  with  the  hot  lilt  jet  lluid 
thickens  the  tree  jc-t  .and  the  w.ill  jet  .iiiti  le.ids  to  a dec'.av  in 
the  m.iximiim  lemper.iture  .inci  veloeilv  with  clist.anee  I roiii  the  lilt 
jet  exit.  The  lie.ited  Wall  jet  will  eventual  Iv  Si-p.ar.lle  I fOIll  t he- 
g.roiind  hei  iiise  ot  hiioy.inl  loic  c s. 


(il  ImpinKi'Hi.'iit  I't  till'  l\i)t  trcf  ji't  ran  appri'c  iah  1 y luMt  t lie  ground 
siirt  ai  l' . 

A solution  of  t lio  oomi>  ross  i b 1 o Koynolds  equations  in  con  j unc  t i c'n  witli 
an  ^ipp  rop  1*  i a t o lurhulonco  mocio!  aiui  thermal  ener^v  ofpiat  ion  coulci  be  used 
to  gain  a f\indamenta!  vindersianding  ol  these  phenomena.  The  results  of 
such  an  analvsis  w<'uld  be  vidoeitv,  pressure,  deiisitv,  and  temperature 
d i St r i but  ions  as  functions  of  the  significant  dimensionless  paramiAers. 
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AITKNDIX  A:  ( ' I'.N'l’KA  I 1)  I I'l'I.K  Al'l’Ki  )X  1 MAT  1 ON'S 

I 'OK  I III'  m'.R  iVAi  1 vi;s 


I'lu  si-i  iMul-i>rdct-. uriir.it  i' , i-i'iii  r.i  1 -cl  i I 1 cri'iui'  .ipiinix  i in;i  t i mis  used  in 
1 tu  ci  i s,  ret  i .'.It  imi  .0  the  <1  i T I e ren  I i .i  I eciu.itimis  are  Kiveii  in  t ti  i s a|)pemiix 
I'in  t i n i t e-(i  i t t e reiii'e  I c > rmii  I ,is  , written  lT>r  ,i  r.eneral  vari.able  , .ire 
expressed  in  terms  at  ['ll  i nt -cii  - 1 lie-ei>mp.i.ss  nol.it  ion  delined  in  Fixture  A- I 
tor  both  interior  .md  bound  irv  point':. 
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Fiqure  A 1 Grid  nol.ition  used  in  the  centr.il  difference  approximation  for  the  derivatives 


lor  ititeriiir  points  .ind  lor  upper-  .md  I owe  r-boimd.i  rv  points: 
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